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ABSTRACT 

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric 
shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness 
of the light curve appears to be compatible both with binary-lens and binary-source models, 
including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. 
The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated 
real-time detection of weak microlensing anomalies with immediate feedback is feasible, effi- 
cient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that 
are not easily resolved from photometric light curves, 3) a modelling approach that finds all 
features of parameter space rather than just the 'favourite model' is required, and 4) the data 
quality is most crucial, where systematics can be confused with real features, in particular 
small higher-order effects such as orbital motion signatures. It moreover becomes apparent 
that events with weak signatures are a silver mine for statistical studies, although not easy to 
exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source 
microlensing events might hide here. 

Key words: gravitational lensing - planetary systems. 



1 INTRODUCTION 



The 'most curious' ef fect of gravitational microlensing jEinsteinl 
ll93d|Pac zvnski"l986^ lets us ex tend our knowledge of planetary 
systems |Mao & Paczviiskj| ll99ll ) to a region of parameter space 
unreachable by other methods and thus populated with intriguing 
suiprises. Microlensing has already impressively demonstrated its 
sensitivity to Super-Earths wit h the detection of a 5 E arth- mass (un- 
certain to a factor two) planet dSeaulieu et al.ll2 006V and it reaches 
down even to about the mass of the Moon l lPacz viiski 1996). 

The transient nature of microlensing events means that rather 
than the characterisation of individual systems, it is the popula- 
tion statistics that will provide the major scientific return of ob- 
servational campaigns. Meaningful statistics will however only 
arise with a controlled experiment, following well-defined fully- 
deterministic and reproducible procedures. In fact, the observed 
sample is a statistical representation of the true population under 
the respective detection efficiency of the experiment. An analysis 
of 13 events with peak magnifications Ao ^ 200 obser ved between 
2005 and 2008 provided the first well-defined sample jGould et al.l 
I2OIOI) . In contrast, the various claims of planetary signatures and 
further potential signatures come with vastly different degrees of 
evide nce and arise fro m different data treatments and applied cri- 
teria jPominikllioiol) as well as observing campaigns following 
different strategies. 

While the selection of highly-magnified peaks is relatively 
easily controllable, and these come with a particiilarlj^ jarge sensi- 
tivity to planetary compa nions to the lens star (G riest & Safizad eh 
1 19981 : 1 Home et al. their rarity poses a fundamental limit 

to planet abundance measurements. Moreover, the finite size of 
the source stars strongly disfavours the immediate peak region for 
planet masses < 10 Af®, where a large magnification results from 
source and lens star being very closely aligned. In contrast, dur- 
ing the wing phases of a microlensing event, planets are more 
easily recognised with larger sources because of an increased sig- 
nal duration, as long as the amplitude exceeds the threshold given 
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by the photometric accuracy (c.f lHaiill2007l) . It is therefore not 
a surprise at all that the two least massive planets found so far 
with unambiguous evidenc e from a well-covered anomaly, namely 
OGLE-2005- BLG-390Lb jSeaulieu et alj|2006h and MOA-2009- 
BLG-266Lb jMuraki etal.ll201ll) . come with an off-peak signature 
at moderate magnification with a larger source star. 

An event duration of about a month and a probability of 
~ 10~^ for an observed G alactic bu lge star to be substantially 
brigh tened at any given time jPaczvrisk i 1991 : Kiraga & Paczyns^ 
1 1994) c alled for a 2-step strat egy of survey and follow-up obser- 
vations i lGould&Loeblll992h . In such a scheme, surveys mon- 
itor > 10^ stars on a daily basis for ongoing microlensing 
events ("Udalski et ai]|l992l : lAIcock et al.|[l99l : iBond et al.ll200ll 
lAfonso et al. 20(3^ whereas roughly hourly sampling of the most 
promising ongoing events with a network of telescopes support- 
ing round-the-clock coverage and photometric accuracy of < 2 % 
allows not only the d etection of planetary signatures, but also 
their characterisation dAlbrow et al.l Il998l: iDominik et J1 12002| : 
[Surgdorf et al. 200 7). While real-time alert syste ms on ongoing 
anomalies l lUdalski et al.ll994l : lAlcock et al.ll99^ . combined with 
the real-time provision of photometric data, paved the way for effi- 
cient target selection by follow-up campaigns, the real-time identi- 
fication of planetary signatures and other deviations from ordinary 
light curves ('anomalies') (IJdalski 2003; Dpminik.et al. 2007) al- 
lowed a transition to a 3-step-approach, where the regular follow-up 
cadence can be relaxed in favour of monitoring more events, and a 
further step of anomaly monitoring at ~ 5-10 min cadence (includ- 
ing target-of-opportunity observations) is added, suitable to extend 
the exploration to planets of Earth mass and below. 

The recent and upcoming increase of the field-of-view of mi- 
crolensing surveys (MOA: 2.2 deg^, OGLE-IV: 1.4 deg^. Wise 
Observatory: 1 deg^, KMTNet: 4 deg^) allows for sampling inter- 
vals as small as 10-15 min. This almost merges the different stages 
with regard to cadence, but the surveys are to choose the expo- 
sure time (determining the photometric accuracy) per field rather 
than per target star. Moreover, they cannot compete w ith the angu- 
lar resolution pos s ible with lucky-imaging cameras (lFriedlll978l : 
iTubbs et al.l kooj : ll0rgensenl I2OO8I : iGrundahl et id] |2009|) . given 
that this technique is incompatible with a wide field-of-view. This 
leaves a most relevant role for ground-based follow-up networks 
in breaking into the regime below Earth mass, in particular with 
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space-based surveys ( iBennett & Rhiel2002h at least about a decade 
away. 

It is rather straightforward to run microlensing surveys in 
a fully-deterministic way, but it is very challenging to achieve 
the same for both the target selection process of follow-up cam- 
paigns and the real-time anomaly identification. The pioneering 
use of robotic telescopes in this field with the Ro boNet cam- 
paigns teurgdorf et"al]|2007l : iTsapras et aLll2009l) led iHome et all 
( |2009|) to devise the first workable target prioritisation algorithm. 
While the OGLE EEWS l lUdalskill2003l) was the first automated 
system to detect potential deviations from ordinary microlensing 
light curves, it flagged such suspicions to humans who would then 
take a decision. In c ontrast, the SIGNALMEN anomaly detector 
jPominik et al.ll2007h was designed just to rely on statistics and 
request further data from telescopes until a decision for or against 
an ongoing anomaly can be taken with sufficient confidence. SIG- 
NALMEN has already demonstrated its power by detecting the first 
sign of fin ite-source effects in MOA-2007-BLG-233/OGLE-2007- 
BLG-302 dChoi et al.|[2oTll) - ahead of any humans - and leading 
to the first alerts sent to observing teams that resulted in crucial 
data being taken on OGLE-2007-BLG-355/MOA-2007-BLG-278 
jHan etalJ |2009|) an d OGLE-2007-BLG-368/MOA-2007-BLG- 
308 dSumi et al.ll201(]|) . the latter involving a planet of ~ 20 Mq. 
SIGNALMEN is now fully embedded into the ARTEMiS (Auto- 
mated Robotic Telesc ope Exoplanet Microlensing Search) software 
system jPominik et a l. 2008b, a) for data modelling and visualisa- 
tion, anomaly detection, and target selection. The 2008 MiNDSTEp 
(Microlensing Network for the Detection of Small Terrestrial Exo- 
planets) campaign directed by ARTEMiS provided a proof of con- 
cept f or fully-deterministic follow-up observations jPominik et al.l 
l2O10h . 

In event OGLE-2008-BLG-510, discussed in detail here, ev- 
idence for an ongoing microlensing anomaly was for the first 
time obtained by an automated feedback loop realised with the 
ARTEMiS system, without any human intervention. This demon- 
strates that robotic or quasi-robotic follow-up campaigns can oper- 
ate efficiently. 

A fundamental difficulty in obtaining planet population statis- 
tics arises from the fact that many microlensing events show weak 
or ambiguous anomalies, sometimes with poor quality data, which 
are left aside because the time investment in their modelling would 
be too high compared to the dubious perspective to draw any def- 
inite conclusions. Indeed such events represent a silver mine for 
statistical studies yet to be designed. Event OGLE-2002-BLG-55 
has already been very rightfully classified as "a possible plane - 
tary event' ' jjaroszvnski & Paczvnskil l2002t iGaudi & Hanll2004l) . 
where ambiguities are mainly the result of sparse data over the sus- 
pected anomaly. Here, we show that OGLE-2008-BLG-510 makes 
another example for ambiguities, which in this case arise due to 
the lack of prominent features of the weak perturbation near the 
peak of the microlensing event. Given that is not a powerful 
discri minator, in pa rticular in the absence of proper noise models 
(e.g. lAnsarilll994l) . it needs a careful analysis of the constraints 
on parameter space posed by the data rather than just a claim of a 
'most favourable' model. In fact, the latter might point to excitingly 
exotic configurations, but it must not be forgotten that maximum- 
likelihood estimates (equalling to minimisation under the as- 
sumption of measurement uncertainties being normally distributed) 
are not guaranteed to be anywhere near the true value. We therefore 
apply a modelling approach that is based on a full classification of 
the finite number of morphologies of microlensing light curves in 



order to make sure that no feature of the intricate parameter space 
is m issed (Bozza et al. , in preparation). 

iDominikI l ll998bh argued that the apparent paucity of mi- 
crolensing e vents reported that involve source rather than lens bi- 
naries (e.g. lGriest&Hu|[l992l) could be the result of an intrin- 
sic lack_o££haractenstic features, but despite a further analysis 
bv lHan&Jeon ^ ^199^) . the puzzle is not solved yet. Moreover, 
while all estimates of the planet abundance from microlensing ob- 
servations indicat e a quite moderate number of mass ive gas giants 
I Sumi et al.|[201ol ; iGould et al.ll201ol ; IPominikllioTTI ; |Cassan et^ 
201 ll) . the small number of reported brown dwarfs (c.f. jPominik 
2010r) . much easier to detect, seems even more striking. As we will 



see, the case of OGLE-2008-BLG-510 appears to be linked to both. 
Maybe the full exploration of parameter space for events with weak 
or without any obvious anomaly features will get us closer to un- 
derstanding this issue which is of primary relevance for deriving 
abundance statistics. 

In Sect. 2, we report the observations of OGLE-2008-BLG- 
510 along with the record of the anomaly detection process and 
the strategic choices made. Sect. 3 details the modelling process, 
whereas the competing physical scenarios are discussed in Sect. 4, 
before we present final conclusions in Sect. 5. 



2 OBSERVATIONS AND ANOMALY DETECTION 

The microlensing event OGLE-2008-BLG-510 at 
RA 18''09'"37:65 and Dec -26°02'26'.'70 (J2000), first dis- 
covered by the Optical Gravitational Lensing Experiment (OGLE), 
was subsequently monitored by several campaigns with telescopes 
at various longitudes (see Table [TJ, and independently detected as 
MOA-2008-BLG-369 by the Microlensing Observations in Astro- 
physics (MOA) team. Table |2] presents a timeline of observations 
and anomaly detection 

When the follow-up observations started, the event m agnifi- 
cation was estimated by SIGNALMEN jPominik et al.ll2007h to be 
A ~ 3.9, which for a baseline magnitude I ~ 19.23 and ab- 
sence of blending means an observed target magnitude / ~ 17.75. 
The event magnification implied an initial sampling interval for the 
MiNPS TEp campaign of r = 60 min according to the MiNPSTEp 
strategy jPominik et aLll2OI0l) . 

The OGLE, MOA, MiNPSTEp, RoboNet-II, and PLANET- 
III groups all had real-time data reduction pipelines running, 
and with efficient data transfer, photometric measurements were 
available for assessment of anomalous behaviour by SIGNALMEN 
shortly after the observations had taken place. While the Micro- 
FUN team also took care of timely provision of their data, we did 
not manage at that time to get a data link with SIGNALMEN in- 
stalled, but since 2009 we enjoy an efficient RSYNC connection^ 
RoboNet data processing unfortunately had to cease due to fire in 
Santa Barbara in early July. After resuming of operations and work- 
ing through the backlog, RoboNet data on OGLE-2008-BLG-510 
were not available before 23 August 2008. 

The SIGNALMEN anomaly detector jPominik et al.ll2007h is 
based on the principle that real-time photometry and flexible 
scheduling allow requesting further data for assessment until a de- 
cision about an ongoing anomaly can be taken with sufficient con- 
fidence. The specific choice of the adopted algorithm comes with 



^ RSYNC is a software application and network protocol for synchronising 
data stored in different locations that keeps data transfer to a minimum by 
efficiently working out differences. — |http : //rsync ■ samba ■ org| 



© 0000 RAS, MNRAS 000, 000-000 



4 V. Bozza et al. 

Table 1. Overview of campaigns that monitored OGLE-2008-BLG-510 and the telescopes used 



Campaign 


Telescope 


Site 


Country 


OGLE 


Warsaw 1.3 m 


Las Campanas Observatory (LCO) 


Chile 


MOA 


MOA l.om 


Mt John University Observatory (MJUO) 


New Zealand 


MiNDSTEp 


Danish 1.54m 


ESO La Silla 


Chile 


MicroFUN 


SMARTS 1.3m 


Cen'o Tololo Inter-American Obsereatory (CTIO) 


Chile 


RoboNet-II 


Faulkes North 2.0m (FTN)* 


Haleakala Observatory 


United States (HI) 




Faulkes North 2.0m (FTS)* 


Siding Spring Observatory (SSO) 


Australia (NSW) 




Liverpool 2.0m (LT) 


Observatorio del Roque de Los Muchachos 


Spain (Canary Islands) 


PLANET-III 


EUzabeth 1.0m 


South African Astronomical Observatory (SAAO) 


South Africa 




Canopus 1.0m 


Canopus Observatory, University of Tasmania 


Australia (TAS) 




Lowell 0.6m 


Perth Observatory 


Australia (WA) 


(ToO) 


IRSF 1.4m 


South African Astronomical Observatory (SAAO) 


South Africa 


(ToO) 


SALT 11m 


South African Astronomical Observatory (SAAO) 


South Africa 



OGLE: Optical Gravitational Lensing Experiment, MOA: Microlensing Observations in Astrophysics, MiNDSTEp: Microlensing Network for the Detection 
of Small Terrestrial Exoplanets, PLANET: Probing Lensing Anomalies NETwork, (ToO): target-of-opportunity observations at further sites not participating 
in regular microlensing observations. 

* The FTN and FTS are part of the Las Cumbres Observatory Global Telescope (LCOGT) Network 



Table 2. Timeline of OGLE-2008-BLG-510 observations and anomaly detection 



Date Time 
28 July 2008 1 5 : 1 2 UT Event OGLE-2008-BLG-5 1 announced by OGLE 

3 August 2008 Event selected by the ARTEMIS system for MiNDSTEp follow-up observations 

4 August 2008 0:50 UT First MiNDSTEp data from the Danish 1.54m at ESO La Silla (Chile) 

1:52 UT First MicroFUN data from the CTIO I.3m (Chile) 

7 August 2008 First RoboNet-II data from the Faulkes North 2.0 (FTN, Hawaii), Faulkes South 2.0m (FTS, Australia), and 

Liverpool 2.0m (LT, Canary Islands) 

8 August 2008 0:39 UT First PLANET-III data from the SAAO 1 .Om (South Africa) 

9 August 2008 1 :30 UT Event announced as MOA-2008-BLG-369 by MOA, following independent detection 

4:58 UT signalmen suspects anomaly based on Danish 1.54m data, acquired at 4:41 UT; 7 further data points were 

taken with this telescope lead to revision of model parameters (peak 8 hours later) 
1 1 :43 UT PLANET-III starts observing the event with the Canopus 1 .Om of the University of Tasmania 
21:01 UT SIGNALMEN suspects anomaly based on SAAO 1 .Om data acquired at 20:28 UT. and with 3 further SAAO 
data points available by 22:00 UT, the model parameters were revised again (peak expected another 8 hours 
later) 

10 August 2008 5:39 UT signalmen triggers check status based on a data point acquired with the Danish 1.54m at 5:36 UT (just 

before dawn) 

6:01 UT Deviation confirmed by further data promptly taken data at 5:47, 5:50, and 5:55 UT, prompting to a stronger 
rise than expected, but SIGNALMEN 1 data point short of calhng an anomaly 

9:09 UT Reassessment triggered by OGLE data point obtained at 3:35 UT, found to be deviating, but final assessment 
unchanged 

9:12 UT The MiNDSTEp and ARTEMiS teams circulate an e-mail to all other microlensing observing teams point- 
ing to an ongoing anomaly 

13:00 UT SIGNALMEN evaluation of 4 data points taken as part of the regular MOA observations between 8:08 UT 
and 10:50 UT led to the automated activation of 'anomaly' status. As a result, further data were taken with 
the telescopes already observing the event, and moreover the IRSF 1.4m and the SALT, both at SAAO 
(South Africa). 

(boldface highlights epochs marked in Fig. |2) 



substantial arbitrariness, where the power for detecting anomalies 
needs to be balanced carefully against the false alert rate. SIGNAL- 
MEN assigns a "status" to each of the events, which is either 'or- 
dinary' (= there is no ongoing anomaly), 'anomaly' (= there is an 
ongoing anomaly), or 'check' (= there may or may not be an on- 
going anomaly). This "status" triggers a respective response: 'or- 
dinary' events are scheduled according to the standard priority al- 
gorithm, 'anomaly' events are alerted upon, initially monitored at 
high cadence, and given manual control for potentially lowering 



the cadence, while for 'check' events further data at high cadence 
are requested urgently until the event either moves to 'anomaly' or 
back to 'ordinary' status. 

SIGNALMEN also adopts strategies to achieve robustness 
against problems with the data reduction and to increase sensitiv- 
ity to small deviations, namely the use of a robust fitting algorithm 
that automatically downweights outliers and its own assessment of 
the scatter of reported data rather than reliance on the reported es- 
timated uncertainties. As a result, SIGNALMEN em on the cautious 
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side by avoiding to trigger anomalous behaviour on data with large 
scatter at the cost of missing pot ential deviation s . The SIGNALMEN 
algorithm has been described bv lDominik et al.l ( l2007h in every de- 
tail. We just note here that suspicion for a deviant data point that 
gives rise to a suspected anomaly is based on fulfilling two cri- 
teria, which asymptotically coincide for normally distributed un- 
certainties: (1) the residual is larger than 95 % of all residuals, (2) 
the residual is larger than twice the median scatter. This implies 
that SIGNALMEN is expected to elevate events to 'check' status for 
about 5 % of the incoming data, but the power of detecting anoma- 
lies outweighs the effort spent on false alerts. In order to allow a 
proper evaluation of the scatter, for each data set, at least 6 data 
points and observations from at least 2 previous nights are required. 
SIGNALMEN moves from 'check' to 'anomaly' mode with a se- 
quence of at least 5 deviant points found. 

On event OGLE-2008-BLG-510, there were two early suspi- 
cions of an anomaly, both on 9 Aug 2008 (see Table |2] Each of 
these lead to a revision of the model parameters with the event ex- 
pected about 8 hours later than estimate earlier, rather than find- 
ing conclusive evidence of an ongoing anomaly. In fact, this be- 
haviour is indicative of the event flattening out its rise earlier than 
expected. Moreover, weak anomalies over longer time-scales look 
marginally compatible with ordinary events at early stages, and fail- 
ure to match expectations can result in a series of deviations that let 
SIGNALMEN trigger 'check' mode, which then leads to a revision of 
the model parameters rather than to a firm detection of an anomaly, 
only to happen once stronger effects become evident. 

On 10 August 2008, SIGNALMEN predicted a magnification 
of A ^ 15.6, which meant a sampling interval of 30 min for 
the MiNDSTEp observations with the Danish 1.54m. Since the 
event was first alerted by OGLE, with a peak magnification Aa = 
4.6 ± 13.30 the signalmen estimate had changed from initially 
Aq ~ 2.3 to ylo ~ 5.3 when the first data with the Danish 
1. 54m were obtained to now Aq ~ 17 (which bears some similar- 
ity with 'model Ic', discussed in the next section). In contrast to a 
maximum-likelihood estimate which tends to overestimate the peak 
magnification, the maximum a posteriori estimate used by SIG- 
NALM EN tends to underestimate it (Albrow 2004; Dom inik et al.l 
l2007h . Just before dawn in Chile, an ongoing anomaly was sus- 
pected again, with further data subsequently leading to firm evi- 
dence. Despite increased airmass likely affecting the measurements 
towards the end of the night, the earlier 'check' triggers were in- 
dicative of a real anomaly being in progress. 

MOA, Canopus 1.0m, and FTN were able to cover the peak 
region, which looked evidently asymmetric. The descent was cov- 
ered by SAAO and then by the Danish 1.54m. Unfortunately, the 
Moon heavily disturbed the observations in the descent, inducing 
large systematic errors that were difficult to correct in the offline 
reduction. Furthermore, we missed three nights of MOA data (12- 
14 August), five nights of OGLE (11-15 August) and we had a two- 
day hole between 12 and 14 August not covered by any telescopes. 
Follow-up observations were resumed after 14 August, including 
PLANET-III observations with the Perth 0.6m (Western Australia) 
from 15 to 20 August, and continued until 30 August. After then, 
only the OGLE and MOA survey telescopes continued to collect 
data. 

We found all data being affected by a scatter larger than what 
might be reasonably expected by the eiTor bars assigned by the re- 
duction software. In fact, following the anomaly around the peak 



The quoted uncertainty refers to a locally lineaiised model. 



of the light curve, all telescopes suffered from inferior data quality 
resulting from increased sky flux contribution due to the proxim- 
ity of the full Moon to the target field as well as contamination by 
a nearby saturated star. The relative faintness of the microlensed 
source (/ = 19.2 according to OGLE database) added to the prob- 
lems. A re-analysis of the images obtained from the FTN, LT, FTS, 
Danish 1. 5 4m, SAAO 1.0m , and Canopus I.Om telescope using 
DanDI^ ( lBramichll2008l) made a crucial difference by remov- 
ing previously present systematics which could have easily been 
mistaken as indications of higher-order effects, and posed a puzzle 
in the interpretation of this event. D anDIA is an im plementation 
of the difference imaging technique jAlard & Luptonlll998i) whose 
specific power arises from modelling the kernel as a discrete pixel 
array, allowing us to properly deal with distorted star profiles be- 
cause it makes no underlying assumptions about the shape of the 
point-spread function. 

In our analysis we have used the data taken at all mentioned 
telescopes except for CTIO, IRSF, SALT, and Perth 0.6m, since 
they are too sparse or too scattered to significantly constrain the fit. 
Moreover, we have neglected all data prior to HJD = 2454500, 
where the light curve is flat and therefore insensitive to the model 
parameters. Furthermore, we have re-binned most of the data taken 
in the nights between 12 and 16 August disturbed by the Moon 
since they were very scattered and redundant. Finally, we have re- 
normalised all the error bars so as to have equalling the number 
of degrees of freedom for the model with lowest y;^, which is re- 
lated to the assumption that it provides a reasonable explanation 
of the observed data. These prepared data sets are shown in Fig. 
[T] along with the best-fitting model that will be presented and dis- 
cussed in detail in Section 4. The peak anomaly is illustrated in 
more detail in Fig.|2] 

Looking at the data and the model light curve, it seems in 
fact that the early triggers on 9 August 2008, spanning the region 
4687.69 ^ HJD - 2450000.0 4688.39 were due to a real 
anomaly, but its weakness together with the limited photometric 
precision and accuracy did not allow us to obtain sufficient evi- 
dence. 



3 MODELLING 

3.1 Microlensing light curves 

Gravitational microlensing events show a transient brightening of 
an observed source star that results from the gravitational bending 
of its light by an intervening object, the 'gravitational lens'. The 
gravitational microlensing effect of a body with mass M is charac- 
terised by the angular Einstein radius 

where G is the universal gravitational constant, c is the vacuum 
speed of light, and 

TTLs = 1 AU (Dl ' - Ds^) (2) 

is the relative parallax of the lens and source stars at distances Dl 
and Ds from the observer, respectively. 

With source and lens star separated by an angle u 9e on the 
sky, the magnification becomes ( lEinsteinlll936h 

DanDIA is built from the DanIDL library of IDL routines available at 
|http : / /www . danidl . co . uk\ 
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UT Date 



22/6 12/7 1/8 21/8 10/9 30/9 20/10 




4620 4640 4660 4680 4700 4720 4740 4760 

HJD-2450000 



Figure 1. Data acquired with several telescopes (colour-coded) on gravitational microlensing event OGLE-2008-BLG-510 (MOA-2008-BLG-369) along with 
the best-titting model light curve and the respective residuals. The region marked by the box is expanded in Fig.|2] 



A{u) 



+ 2 



u ^u^ + 4 



(3) 



If we assume a uniform relative proper motion \i between lens and 
source star, the separation parameter u becomes 



u{t;to,uo,tE) = ^ 



(4) 



where Ie = 6'e//i is the event time-scale, and the closest angular 
approach uq is realised at time to- 

With Fs being the flux of the observed target star, and Fb the 
background flux, the total observed flux becomes 

F{t) = FsA{t) +Fb= Fbaso ^l^l^'^ = FbaBcAobs(t) (5) 

with the baseline flux Fbasc = Fs + Fb and the blend ratio g = 
Fb/Fs, where 



+ 9 
1+5 



(6) 



is the observed magnification. 

Because of A{u) monotonically increasing as u — s> 0, or- 
dinary microlensing light curves (due to single point-like source 
and lens stars) are symmetric with respect to the peak at to, where 
the closest angular approach between lens and source u(to) = uo 
is realised, and fully characterised by (fo, ito,fE) and the set of 
(Fs,Fb) or (Fbasc, (?) for each observing site and photometric 
passband. Best-fitting (Fs, Fb) follow analytically from linear re- 
gression, whereas the observed flux is non-linear in all other pa- 
rameters. 



3.2 Anomaly feature assessment and parameter search 

Apparently, the only evident feature pointing to an anomaly in 
OGLE-2008-BLG-510 is the asymmetric shape of the peak (see 
Fig |2]l. Such weak effects can be explained by a finite exten- 
sion of the central point caustic of a single isolated lens star due 
to binarity (which includes the presence of an orbiting planet). 
Moreover, the absence of strong effects excludes the source star 
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Figure 2. Data and model light curve for OGLE-2008-BLG-510 around the asymmetric peak. For comparison, we also show best-fitting model light curves 
for a single lens stai' using all data (short dashes) or data before triggering on the anomaly on 10 August 2008 (long dashes). Moreover, we have indicated the 
most relevant stages in the real-time identification of the anomaly. 



from hitting or passing over the caustic. This straightforwardly 
restricts parameter space, and one could readily identify a lim- 
ited number of viable configurations with regard to the topolo- 
gies of the caustic and the source trajectory and exclude all oth- 
ers. The use of an 'event library' where the most important fea- 
tures are stored had already been suggested by Mao & Di Stefano 
(Il995h . while generic features were the starting points for the ex- 
ploration of para meter space in early efforts of modelling mi - 
crolensing events ( lDomimir& Hirshfeldl [l996l ; iDominikI 1 1 999ah . 
Moreover, the identification of features for caustic-passage events 
has been pr oven poweiful in efficient searche s for all viable con - 
figurations jAlbrow et all 1 19991; ICassaiil 120081 ; iKains et al.ll2009h . 
More recent work built upon the univ ersal topologies of binary- 



curves (Bozza 



lens sy stems i Erdl & Schneider 



2001 



'; lNight 



et al. 



19931 



l200i 



in order to classify light 
Based on the earlier work 



by iBozza ( 20Q1 ), we adopted an automated approach (Bozza et 
al., in preparation) that starts from 76 different initial conditions 
covering all possible caustic crossings and cusp approaches in all 



caustic top ologies occurring in binary lensing (close, intermediate 
and wide) ferdl & Schneideilll993l; lDominiki[l999bl) . From these 
initial conditions, we have run a Levenberg-Marquardt algorithm 
for downhill fitting, and then we have refined the "x^ minima by 
Markov chains. 



The roundish shape of the peak however does not allow us 
to immediately dismiss the alternative hypothesis that the source 
rather than the lens is a binary system. Binary-source light curves 
are si mply the superposition of ordinary light curves dOriest & Hij 
h992h . leading to significant zoo of morphologies, which is how- 
ever less diverse than that of binary lenses. Binary-lens systems can 
be uniquely identified from slope discontinuities and the sharp fea- 
tures that are associated with caustics, while smooth, weakly per- 
turbed light curves may be ambiguous, and even potential pl anetary 
signa tures might be mimicked by binary-source systems dOaudil 
h998h . 
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V l+q ) \\+q ) 

Figure 3. The geometry of a binary gravitational lens being approached 
by a single source star and the related parameters. The separation be- 
tween the primary and the secondary is given by d, while the mass ratio is 
q = A/2 /Ml. Positions on the sky aiise from multiplying the dimension- 
less coordinates with the angular Einstein radius . so that = ''e /a* 
given an event time-scale, where ^ denotes the absolute value of the relative 
proper motion between source and lens star 

3.3 Binary-lens models 

In addition to its total mass M = Mi + M2, a binary lens is fully 
characterised by the mass ratio q = A/2/A/1 of its constituents, 
and, if one neglects the orbital motion, by their separation. The lat- 
ter can be described by the dimensionless parameter rf, where d 6^ 
is the angle on the sky between the primary and the secondary as 
seen from the observer. In contrast to a single lens, the microlens- 
ing light curve depends on the orientation of the source trajectory, 
where we measure the angle 6 from the axis pointing from M2 to 
Ml. As the reference point for the closest angular approach be- 
tween lens and source, characterised by to and ito, we choose the 
centre of mass. This means that the source trajectory relative to the 
lens is described by 

, , / — sin 6 \ t — to f cos 6 \ 

\ cos 6 J tE \ sm9 J ' 

while the primary of mass Mi is at the angular coordinate [dq/{l + 
q) , 0] Oe and the secondary of mass M2 at the angular coordinate 
[-d/(l + g),0] 6»E (see Fig.[3j. 

Strong differential magnifications also result in effects of the 
finite size of the source star on the light curve, quantified by the 
dimensionless parameter p*, where 6e is the angular source ra- 
dius. We initially approximate the star as uniformly bright. 

For the evaluation of the magnification for given mod el param- 
eters, we have adopted a contour integration algo rithm dPominikl 
I1995I : iGould & Gaucherelll997l : lDominiklll998al) improved with 
parabolic correction, optimal sampling and accurate error esti- 
mates, as described in detail bv lBozzal OOlOh . 

Our morphology classification approach leads to three viable 
configurations where the source trajectory approaches the caustic 
near the primary with different orientation angles, grazing one of 
the four cusps before having passed a neighbouring cusp at larger 
distance. We r ecover the well-known ainbiguity between close and 
wide binaries tGriest & S afizadeh 1998; Dominik 1999b): all con- 
figurations come in two flavours. Fig. |4] illustrates these config- 
urations, labelled T, '2', and '3', for the close-binary topology, 
whereas the wide-binary case is analogous. Given the symmetry of 
the binary-lens system with respect to the binary axis, there are 3 
different cusps, one off-axis and two on-axis. For the off-axis cusp, 
there are two different neighbouring cusps, distinguishing models 



2 and 3. In contrast, the neighbouring cusps to an on-axis cusp is 
the identical off-axis cusp, so that model 1 is not doubled up. In 
principle, there could be a solution near the approach of the other 
on-axis cusp, but this turned out not to be viable. 

We end up with six candidate models that we label by Ic, 2c, 
3c, Iw, 2w, 3w. The "c" con'esponds to the close-binary topology 
and the "w" con'esponds to the dual wide-binary topology. The val- 
ues of the parameters of these models, with their uncertainties are 
shown in Table [S] Model Ic comes with the smallest (set to 
the number of degrees of freedom by rescaling the photometric un- 
certainties), with model Iw closely following with just Ax^ ~ 1- 
Models 2 and 3 come with A^^ ~ 20, with model 2w being sin- 
gled out by its wide parameter ranges. Models 1 prefer an OGLE 
blend ratio close to zero, whereas models 2 prefer a larger back- 
ground, but are still compatible with zero blending. Models 3 come 
with a negative blend ratio. If one imposes a non-negative blend 
ratio, models I and 2 change rather little (see Table |4), but we did 
not find corresponding minima for the configurations of models 3, 
which rather tend towards models 2. Mass ratios q roughly span 
the range from 0. 1 to 1 . Given that the source passes too far from 
the caustic to provide significant finite-source effects, all models 
return only an upper limit on the source size parameter p* . We will 
compare all models in detail in Sect. l3.7l 

3.4 Binary-source models 

As anticipated, the binary source interpretation might be able to 
explain the anomaly in OGLE-2008-BLG-5I0 as efficiently as the 
binary lens. If we neglect the orbital motion, the microlensing light 
curve of a "static" binary source with a single lens is just the super- 
position of two Paczyhski curves with the same ts, so that 

A{t) = {1 - Lj) A[u{t-ti,Ul,tE)] + CO A[u{t-t2,U2,tE)] , (8) 

where i^ = Fs.2/(Fs.i+Js,2)is the flux offset ratio for the source 
fluxes Fs.i and Fs.2, and we just need to distinguish two pairs 
{ui,ti) and {u2,t2) that characterise the closest angular approach 

to each of the constituents^ 

As pointed out bv lDominik & Hirshfel j ( Il996l Appendix C), 
there is a two-fold ambiguity for the angular separation At; be- 
tween the source stars, given that the photometric light curve does 
not tell us whether they are on the same side ('cis' configuration) 
or on opposite sides ('trans' configuration) of the source-lens tra- 
jectory, where 



A77 = ^E y (^^^^j +{u2±mf. (9) 

The search in the parameter space is much simpler, since there 
is only one way of superposing two Paczyiiski curves so as to obtain 
the asymmetric peak of OGLE-2008-BLG-510. We found that a 
small negative blend ratio is preferred, where the parameters shift 
only very little if a non-negative blend ratio is enforced. The best- 
fitting model with this constraint is given in Table|5] 

3.5 Higher order effects 

Beyond the basic static binary-lens and binary-source models pre- 
sented above, we looked for potential signatures of three possible 
higher order effects: annual parallax, lens orbital motion and source 
limb darkening. 

As the Earth revolves around the Sun, the trajectory of the 
source relative to the lens system effectively becomes curved. The 
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Table 3. The six static binary lens models explaining the peak anomaly. dB^, is the binary angulai' separation, and q = M2/A/1 the mass ratio. The closest 
angular approach uq 9^ between the source and the centre of mass of the binary lens occurs at time to . The angle 9 measures the orientation of the source 
trajectory from the binary axis, and is the event time-scale. Finally, p^, 9^ is the angular radius of the source star, and g the blend ratio, i.e. the quotient 
of background and source flux, is in units of days, to is HJD— 2450000, 9 is in radians, and all other parameters ai'e dimensionless. We have rescaled the 
photometric uncertainties so that the of model Ic matches the number of degrees of freedom (568). The quoted parameter intervals correspond to levels 
that include 68 % of the Markov Chain realisations; for p* we only give an upper limit at 68 % confidence level. While a non-negative OGLE blend ratio is 
compatible with models 1 and 2, models 3 need adjustment. 
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Table 4. Viable static binary lens models with non-negative OGLE blend ratio imposed. While there is little effect on models 1 and 2, models 3 have dropped 
out with no corresponding minima found. 



the source trajectory characterise the relative direction of the eclip- 
tic. A typical signature of the annual parallax effect in microlensing 
light curves is an asymmetry between the rise and the descent. 

Starting from the static solutions of Tab. |4] we have run 
Markov chains including the two parallax parameters, where our 
results are summarised in Tab.|6] For models 1, reduces by less 
than one with the two additional parameters, which shows the in- 
significance of parallax. There is a larger Ax^ for models 2, namely 
7.3 and 2.5, respectively, but we need to consider that improve- 
ments at such level can easily be driven by data systematics (not fol- 
lowing the assumption of uncorrelated normally-distributed mea- 
surements). We therefore only use an upper limit on the parallax 
for constraining the physical nature of the lens system. Similarly, 
we do not find any significant parallax signal with the binary-source 
model (see Table|5j and a similarly weak limit. 

The orbital motion of the lens might particularly affect our 
candidate models in close-binary topology, where a relatively short 
orbital period compared to the event time-scale Ib is possible and 
worth being checked for. We have therefore performed an extensive 
search for candidate orbital configurations, including the three ve- 
locity components of the secondary lens in the centre of mass frame 
as additional parameters. Due to the lack of evident signatures of 
orbital motion, we have just considered circular orbits, which are 
completely determined by the three parameters just introduced. In- 
deed, we find particular solutions with a sizeable improvement in 
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Table 5. Binary-source model without and with parallax, where we con- 
strain the OGLE blend ratio to be non-negative. The times ti and ta are in 
HJD — 2 450 000, while all other parameters are dimensionless. 



curvature depends on the orientation of the source trajectory rela- 
tive t o the ecliptic , which can be expressed by a parallax vector tte 
(e.g. lGouldl2004) with the absolute value 

(10) 

where the components {'ite,\\, ''^e.±) parallel and perpendicular to 
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Table 6. Models lc+, lw+, 2c+, and 2w+ including the annual parallax effect. The parameters are described in Tab. |3] with the addition of the parallax 
parameters ttj; y and TrE,^. For the total parallax tte we give the upper limit at 68% confidence level. 



the with respect to all static models. For example, model Ic gets 
down to = 539. However, a closer look at the candidate so- 
lutions thus found shows that they are actually fitting the evident 
systematics in the data taken during the descent of the microlens- 
ing light curve, when the Moon was close to the source star. In all 
these data, a weak overnight trend is present, whose steepness de- 
pends on the particular image reduction pipeline employed. Again, 
we see that systematics in the data prevent us from drawing firm 
conclusions just from improvements in x^^ not knowing whether 
the data related to it show a genuine signal. Consequently, we will 
withdraw the hypothesis that orbital motion can be detected, and 
we will not use any orbital motion information in the interpretation 
of the event. 

We have also looked into source orbital motion, and not that 
surprisingly, we also find that the fit is driven by systematic trends 
of the data over the night. Therefore, we will not consider the im- 
provement in the x^ down to 544 as due to intrinsic physical effects. 

Finally, we have tested several source brightness profiles to 
account for limb darkening for all models, but the difference with 
the uniform brightness models is tiny, orders of magnitude smaller 
than the photometric uncertainty. This is consistent with the fact 
that in all models the source passes relatively far from the caustic, 
leaving no hope to study the details of the source structure. 



3.6 Limitations due to systematics in the data 

The photometric analysis for this event posed several challenges. 
The source star was very faint (J ~ 19.23) and heavily blended 
with a brighter companion. A saturated star was nearby, which was 
necessarily masked by our photometric pipeline, but which in turn 
placed constrains on the size of the PSF and the kernel that could be 
used during the image subtraction stage, particularly for the images 
that had high seeing values. In a few specific cases, using a larger 
fit radius would mean that the masked area around the saturated 
star fell within the fit radius of the source s tar, thereby reducing the 
number of pixels used in the photometry ( ISramich et al.||201 ll) . A 
new version that is cuixently under development discounts a bigger 
region of pixels around the saturated star from the kernel solution, 
since the kernel solution is most sensitive to contaminated pixels. 

Shortly after the ongoing anomaly had been identified, the 
moon was full and close to the target field of observation result- 
ing in high sky background counts and a strong background gra- 
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Figure 4. Illustration of the viable caustic and source trajectory configura- 
tions, showing how the trajectory approaches a cusp and passes with respect 
to the caustic. Models 1 and 2 differ in the exchanged roles of the on-axis 
and off-axis cusp, whereas models 2 and 3 exchange the on-axis cusp. Ge- 
ometry and symmetry would allow for a 4th class of models, where the 
closest cusp approach is to the on-axis cusp on the side opposite the com- 
panion, but such were found not to be viable. 
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Figure 5. Observable magnitude shift Am = 2.5 lg{[^(f) + + as a function of time for all suggested models as compared to model Ic. The OGLE 
blend ratio has been used as reference. 



dient. While the pipeline can account for the latter, the former has 
an impact on the photometric accuracy that can be achieved. In 
fact, the moon started to have a strong effect from 1 1 August 2008 
(illumination fraction 80%) and ~7-10° from the target field, on 
to 12 August 2008 (illumination fraction 86%) and ~ 3-7° from 
the target field. It was full on 16 August 2008. So all the affected 
data were taken after the anomaly at the peak. Normally it is not 
advisable to observe if the moon is bright and closer than 15° to 
the target. Obtaining a full characterisation of the impact of Moon 
pollution on the photometry is not that straightforward to assess 
with difference imaging, because there are many other factors that 
affect the kernel solution as well, and the different contributions 
are not easy to isolate. We have however optimised our photometry 
to minimise the impact of Moon pollution and accounted for the 
uncertainty during our modelling runs. Moreover, DanDIA takes 
care of systematic effects by producing a value for the star fits 
and reflecting this in the reported size of the photometric error bars. 



3.7 Comparison of models 

For the peak region, the differences between the models remain be- 
low 2.5 % (see Fig. O, which makes these difficult to probe with 
the apparent scatter of the data. Moreover, only if it is ensured that 
systematic effects are consistently substantially below this level, 
will a meaningful discrimination between the competing models 
be possible. In our opinion, this poses the most substantial limit to 
our analysis, and we therefore abstain from drawing very definite 
conclusions about the physical nature of the OGLE-2008-BLG-510 
microlens. Fig. |5] also shows that the differences between the re- 
spective pairs of close- and wide-binary models are even much 
smaller, less than 0.4 to 0.7 % over the peak. It can also be seen 
that models are made to coincide better in the peak region where 
more data have been acquired, as compared to the wings where 
larger differences of 3-6 % are tolerated. 

As Fig.[6]shows, the differences between the models with re- 
spect to the residuals with observed data appear to be rather sub- 
tle. Additional freedom that causes further narrowing arises from 
adopting a blend ratio for each of the sites and each of the models, 
allowing for different relations between the blend ratios between 
the models. A successful fit should be characterised by data falling 
randomly to both sides of the model light curve. On this aspect, 
models 2 appear better balanced between all data sets than mod- 
els 1, in particular with respect to the Canopus 1.0m and the MOA 
data, but it introduces a trend in some SAAO data over the peak 



that does not show with other models. We also see that models 3 
are more closely related to models 2 than to models 1, with rather 
similar behaviour over the peak, while the differences in the wing 
regions are related to the different blend ratios. In fact, both mod- 
els 2 and 3 have the second and closer cusp approach to the off-axis 
cusp, whereas the roles of the on-axis and off-axis cusps are flipped 
for models 1 . 

It is also very instructive to look at Ax^ between the mod- 
els as a function of time as the data were acquired (see Fig. |7}, 
which obviously depends on the sampling. While models 3 as well 
as the binary-source model rather gradually lose out to models 1 
for 4683 ^ HJD - 2450000 ^ 4691, models 2 perform much 
worse than all other models for the rather short period 4688.0 ^5 
HJD-2450000 ^ 4689.5, which happens just to coincide with the 
skewed peak. However, for 4686.0 ^ HJD - 2450000 4688.0 
and 4689.5 ^ HJD - 2450000 s: 4690.5, models 2 do a better 
job than models 1. It is in the nature of minimisation that more 
weight is given to regions with a denser coverage by data. As a 
consequence, relative between competing regions of parameter 
space depend on the specific sampling. In fact, the region 4689.5 ^ 
HJD — 2450000 ^ 4690.5 which favours model 2 has a sparser 
coverage than preceding nights. Moreover, we find model 3c out- 
performing all other models for 4645 =^ HJD - 2450000 ^ 4682, 
but this is given little weight. The weight arising by the sampling 
is a particularly relevant issue given that for non-linear models, the 
leverage to a specific parameter strongly depends on the time ob- 
servations are taken. We immediately see that there is an extreme 
danger of conclusions being driven by systematics in a single data 
set during a single night, potentially overruling all that we learn 
from other data. 



4 PHYSICAL INTERPRETATION 

An ordinary microlensing event is determined by the fluxes of the 
source star Fs, the blend Fb, the mass of the lens star M, the source 
and lens distances Ds and Dl, as well as the relative proper motion 
/i between lens and source star (or the corresponding effective lens 
velocity v — Di, /i), and finally the angular source-lens impact 
6o and its corresponding epoch to. However, the photometric light 
curve is already fully described by a smaller number of parameters. 
Apart from Fs, and Fb, it is completely characterised by to,uo = 
60/ 6e and = ^e/m, where the angular Einstein radius 6e, as 
given by Eq. Q absorbs M, Di,, and Ds- All the physics of Ds, 
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Figure 6. Residuals in the peaJc region of OGLE-2008-BLG-5 10 for all proposed models. 
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Figure 7. Ax^ as a function of time for all suggested models relative to 
model Ic, along with the cumulative number of data points up to that time. 



Dl, M, and v gets combined into the single model parameter t-^, 
and information about the detailed physical nature gets lost. 

Higher-order effects can play an important role in recovering 
the information. Additional relations between the physical prop- 
erties can be established if the light curve depends on further pa- 
rameters that are related to the angular Einstein radius 6'e or the 
relative lens-source parallax ttls- In particular, finite-source effects 
with the parameter = 6i,/6^ or annual parallax with the mi- 
crolensing parallax tte = ttls/^e allow solving for Dl, M, and 
w, once Ds and the angular size 6', of the source star can be es- 
tablished. Moreover, as in the case of OGLE-2008-BLG-510 here, 
even the absence of detectable effects can provide valuable con- 
straints to parameter space. In fact, as discussed in the previous 
section, the anomaly of OGLE-2008-BLG-510 is almost indiffer- 
ent to the source radius. Moreover, no parallax signal is clearly de- 
tected in either model. Finally, our attempts to find solutions with 
orbital motion just bumped into systematic overnight trends in the 
data. 

But even if the physical properties of the lens system cannot be 
determined directly, Bayes' theorem provides a means of deriving 
their probabilty density. In fact, with the physical properties -0 and 
the model parameters p, one finds a probability density P-,p(xf>\p) 
of -0 given p as 



n(V'lp) 



£(p|^)P^(V>) 



(II) 



where Pif,{tl!) is the prior probability density of the physical prop- 
erties tp, and C{p\\l>) is the likelihood for the parameters p to arise 
from the properties xp. 

While prior probability densities (tp) for the physical prop- 
erties are straightforwardly given by a kinematic model of the 
Milky Way and mass function of the various stellar populations, the 
likelihood £(p|'!/j) is propor tional to the dif ferential microlensing 
rate d'T/(dpi . . . dpk) (c.f. lDominikll2006h . which is to be evalu- 
ated using the constraining relations between the model parameters 



p and the physical properties -0. We determine P^(ip) follo wing 
the lines of Domiiiik ( 2006). While lCalchi Novati et ailoOOSl) have 
recently discusse d Galactic models in so me detail, we basically fol- 
low the choice o f iGrenacher eTZlJigggl) . In particular, we consider 
candidate lenses in an exponential disk or in a bar-shape d bulge . 
The respective mass functions are t aken from IChabrieij (|2003|) . 
The only difference with respect to IPominikl ( |2006|) is that we 
consider an anisotro pic velocity dispersion for stars in the bulge 
l lHan&Goul3l 19951) , with cr = (116, 90, 79) km/s along the three 
bar axes X' ,Y' , Z' defined therein. Actually, this difference does 
not have any major effects on the final result. The denominator in 
Eq. dl It just reflects that the probability density is properly nor- 
malised, i.e. J P-^ ('/''Ip) d0' = 1, which can be achieved trivially. 

In addition to the information contained in the event time- 
scale fE, we also take into account the constraints arising from the 
upper limits on the source size, the parallax, and the luminosity of 
the lens star, as given by the blend ratio g. 

The angular source ra dius can be measured using the tech- 
nique of lYooetal.1 l l2004h . The dereddened color and magnitude 
of the source are measured from an instrumental colour-magnitude 
diagram as [{V — 7),7]o,s = (0.63, 17.78), where we assumed 
that the clump centro id is at [(V - J)o, Mj]ciump = (1.06,-0.10) 
( lBensbvetaUl20Ill ; D. Nataf, private communication), that the 
Galactocentric distance is _Ro = 8.0 kpc and that the clump cen- 
troid lies 0.2 mag closer to us than the Galactic centre, since the 
field is at Galactic latitude I — +5.2°. Therefore, we also set 
Ds ~ 7.3 kpc. W e then convert from V/J to V/K using the colour- 
colour relations of lBessell & Bretll(ll988h. Finallv. we use th e V/K 
colour/surface-brightness relation of Kervella et al.l ( |2004|) to ob- 
tain 9i, = 0.80 \ias. 

Rather than just basing our analysis on the best-fitting values, 
we also take into account the model parameter uncertainties. From 
our Markov-Chain Monte Carlo runs, we find that ts, p*, tte and g 
(OGLE) are basically uncorrelated. Their combined probability can 
therefore fairly be approximated by the product of the individual 
probability densities 

PtE,P*,nE.9{iE,p*,7VE,g) = PtE{tE,)Pp^{p*)P-^E{'^E)Pg{g) .(12) 

We approximate each of these distributions by a Gaussian 

1 (^-s:f 

G(x,x,a) — :=e 



(13) 



where 



-PtE(*E) = G{tE,tE,CrtE) 

Pp,(p.) = 2e(/a.)G(/a*,0,ap,; 
P^j,(7rE) = 2e(7rE)G(7rE,0,o-, 
P,(,g) = 2e(5)G(,g,0,a,), 



(14) 



with the dispersions of these Gaussians being identified with the 
68% confidence limits listed in Table|6] 

We find a probability density in _Dl, M, and v. By integrating 
over V, and changing from M to logM, we obtain a probability 
d ensity in the (Dl , loK M) plane. In this integration, as explained 
in lPominiS ( |2006|) . the dispersion in the Einstein time is practically 
irrelevant and we can approximate the Ie distribution by a delta 
function. The arising posterior probabilities are illustrated in Fig. 
[8] separately for disk and bulge lenses for models 1 and 2. We can 
appreciate the small difference between the close- and wide-binary 
models (whose contours are dashed in Fig. [8). The wide topology 
slightly favours a smaller distance and a higher mass for the lens 
system. Higher masses are cut off by the blending constraint, as the 
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Table 7. Expectation values of lens distance Dl and lens mass M for the 
different binary-lens models, and assessment of probability for the lens to 
reside in the Galactic disk or bulge. 



distance Ds becomes Ds 6'e = 2.35 AU. For a face-on circular 
orbit, we therefore find the orbital radius simply as ro ~ DsArj, 
where Ary as given by Eq. l|9j denotes the angular separation be- 
tween the binary-source constituents. With the total mass of the 
binary-source system m = mi + m2 ~ 1.83 Mq, we then ob- 
tain orbital periods of Po.cis ~ 16 days or Po. trans ~ 29 days 
for the cis- or trans-configurations, respectively. For an orbital in- 
clination i, the orbital radius increases as r — ro/(cosi), so that 
the orbital periods increase as P = Po/(cosi)^''^. For orbits with 
an eccentricity e, the observed separation r is limited to the range 
a(l — e) ^ r ^ a(l + e), where a denotes the semi-major axis. 
While r = a for a circular orbit, one finds in particular a ^ r/2. 
Therefore, in the general case, orbital periods still have to fulfil 
P ^ Po/i'i cosi)^^^. Given both the small colour difference and 
a plausibly long orbital period (as compared to the event time-scale 
of tE ~ 22 days), there is no apparent inconsistency within the 
proposed binary-source model. 



lens would become too bright to be compatible with the absence of 
blending. In fact, this blending constraint readily dismisses mod- 
els 3. The source radius constraint slightly cuts small lens masses 
as these would yield a too small Einstein radius. Finally, only the 
tail of the disk distribution at small Dl is affected by the parallax 
constraint. 

The relative weight of disk and bulge for each microlensing 
model is obtained by integrating these distributions over the whole 
(DL,logM) plane. Table |7] shows the probabilities for each mi- 
crolensing model and each hypothesis for the lens population. We 
also show the average lens distance and mass in each case. We find 
that the cut on higher masses due to the absence of blending sig- 
nificantly penalises the bulge, which is known to have a higher mi- 
crolensing rate thanks to massive lenses. This fact witnesses how 
the information coming from the study of our specific microlens- 
ing event effectively selects the viable lens models. 

Within the parameter uncertainties of the binary-lens models, 
there is some significant overlap in the possible physical nature of 
the lens system (with model 2w being compatible with a very wide 
range). We find that an M dwarf orbited by a brown dwarf appears 
to be a favoured inteipretation, but we cannot fully exclude a more 
massive secondary. Regardless of the specific model, a location in 
the Galactic disk appears to be as plausible as in the Galactic bulge. 
The microlensing light curve does not definitely clarify whether 
the system is in a wide or in a close configuration because of the 
well-known degeneracy of the central caustic. It is interesting to 
note that if we use the average values for the total mass and the 
distance we can estimate that the minimum value for the orbital 
period in model Ic (Iw) is 296 days (51 years) with a minimum 
orbital radius of 0.52 (8.25) AU. Even in the closer configuration, 
detecting any orbital motion signal in this event would have been 
quite difficult. This also reinforces our rejection of all models with 
very short orbital periods found during the modelling phase. 

For the binary-source model, we find with the I band magni- 
tude and the luminosity offset ratio u, the masses of the two com- 
ponents as mi = 1.03 Mq and m2 = 0.80 Mq. The V ~ I 
index of this system is 0.1 mag larger than the V — I index of 
a single source that would yield all the observed flux. Moreover, 
by adopting values for the lens mass and distance in the middle 
of the two averages for the bulge or disk hypothesis (see Table |7]l, 
the physical size of the angular Einstein radius 9b at the source 



5 CONCLUSIONS 

So far, most of the results arising from gravitational microlensing 
campaigns have been based on the detailed discussion of individual 
events that show prominent characteristic features that allow unam- 
biguous conclusions on the underlying physical nature of the lens 
system that caused the event. However, with its sensitivity to mass 
rather than light, gravitational microlensing is particularly suited 
to determine population statistics of objects that are undetectable 
by any other means. This not only includes planets in regions inac- 
cessible to other detection techniques, but also brown dwarfs, low- 
mass stellar binaries, and stellar-mass black holes. Meaningful pop- 
ulation statistics will only arise from controlled experiments with 
well-defined deterministic procedures, and similarly the data anal- 
ysis needs to adopt a homogeneous scheme. Given that the infor- 
mation arises from the constraints posed by the data on parameter 
space, i.e. what can be excluded rather than what can be detected, 
events with weak and potentially ambiguous features need to be 
dealt with appropriately. 

The event OGLE-2008-BLG-510, which was subject to a de- 
tailed analysis in this paper, is a rather typical representative of 
the large class of events with evident weak features. Our light 
curve morphology classification approach revealed that the appar- 
ent asymmetric shape of the peak appears to be compatible with 
a range of binary-lens models or a binary-source model. A par- 
ticular challenge in the interpretation of the acquired data is deal- 
ing with systematic effects that can lead to the mis-estimation of 
measurement uncertainties, non-Gaussianity, and correlations. This 
in turn strongly limits the validity of the application of statistical 
techniques that assume uncorrelated and normally (Gaussian) dis- 
tributed data. We moreover find that the sampling of microlensing 
light curve affects the preference of regions of parameter space, and 
that conclusions therefore are at risk to be driven by systematics in 
the data rather than real effects. In the specific case of OGLE-2008- 
BLG-510, all suggested models coincide within 2.5 % in the peak 
region (where the sampling is dense), and in order to discriminate, 
not only does the scatter of the data need to be small, but more im- 
portantly, systematics need to be understood or properly modelled 
to a level substantially below the model differences. For OGLE- 
2008-BLG-510, we found that systematic trends in some datasets 
simulate a fake scale for orbital motion, which (if present) is hidden 
below systematics and thus impossible to determine. 

Despite the fact that the viable binary-lens models differ in the 
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Figure 8. Contours of the probability density containing 95%, 68%, 38%, 10% probability in the mass-distance plane for the lens in binary-lens models 1 
(top), binary-lens models 2 (middle), and the binary-source model (bottom). The panels on the left treat the case of a disk lens, while those on the right consider 
a bulge lens. The dashed lines are the contours for models Iw and 2w, respectively. The filled (empty) circles are the modes of the probability density for the 
close- (wide-)binary models. 



geometry between the binary-lens system and the source trajectory, 
there is some overlap in the resulting physical nature of the binary- 
lens system, taking into account the model parameter uncertainties. 
In fact, we find the lens system consisting of an M dwarf orbited 
by a brown dwarf a favourable interpretation, albeit that we cannot 
exclude a different nature. 

The large number of microlensing events with weak features 
are currently waiting for a comprehensive statistical analysis, and 
lots of interesting results are likely to be buried there. It however 
requires a very careful analysis of the noise effects in the data as 
well as modelling techniques that map the whole parameter space 
rather than just finding a single optimum, which is just a point esti- 
mate, and moreover biased if minimisation is usedQ Currently, 

^ This already holds if all data ai'e uncorrelated and unifonnly distributed. 
Any deviation from these assumptions makes it even worse. 



we neither understand the paucity of reported events due to a binaiy 
source nor the impact of brown dwarfs on microlensing events, and 
we might even need to look into weak signatures and their statistics 
in order to arrive at a consistent picture. 



The detection of the weak microlensing anomaly in OGLE- 
2008-BLG-510 also demonstrated the power and feasibility of au- 
tomated anomaly detection by means of immediate feedback, i.e. 
requesting further data from a telescope which is automatically re- 
duced and photometric measurements being made available within 
minutes. This is an important step in the efficient operation of a 
fully-deterministic campaign that gets around the involvement of 
humans in the decision chain and offers the possibility to derive 
meaningful population statistics by means of simulations. 



© 0000 RAS, MNRAS 000, 000-000 



16 V. Bozza et al. 



ACKNOWLEDGMENTS 

The Danish 1.54m telescope is operated based on a grant from the 
Danish Natural Science Foundation (FNU). The "Dark Cosmol- 
ogy Centre" is funded by the Danish National Research Founda- 
tion. Work by C. Han was supported by a grant of National Re- 
search Foundation of Korea (2009-0081561). L. M. acknowledges 
support for this work by research funds of the International Insti- 
tute for Advanced Scientific Studies. Work by AG was supported 
by NSF grant AST-0757888. Work by BSC, AG, and RWP was 
supported by NASA grant NNX08AF40G. Work by SD was per- 
formed under contract with the California Institute of Technology 
(Caltech) funded by NASA through the Sagan Fellowship Program. 
The MOA team acknowledges support by grants JSPS20340052, 
JSPS20740I04 and MEXT19015005. Some of the observations re- 
ported in this paper were obtained with the Southern African Large 
Telescope (SALT). LM acknowledges support for this work by re- 
search funds of the International Institute for Advanced Scientific 
Studies. MH acknowledges support by the German Research Foun- 
dation (DFG). DR (boursier FRIA) and JSurdej acknowledge sup- 
port from the Communaute franfaise de Belgique - Actions de 
recherche concertees - Academie universitaire Wallonie-Europe. 
The OGLE project has received funding from the European Re- 
search Council under the European Community's Seventh Frame- 
work Programme (FP7/2007-2013) / ERC grant agreement no. 
246678. MD, YT, DMB, CL, MH, RAS, KH, and CS are thank- 
ful to Qatar National Research Fund (QNRF), member of Qatar 
Foundation, for support by grant NPRP 09-476-1-078. 



REFERENCES 

Afonso C, et al., 2003, A&A, 404, 145 
Alard C, Lupton R. H., 1998, ApJ, 503, 325 
Albrow M., et al., 1998, ApJ, 509, 687 
Albrow M. D., 2004, ApJ, 607, 821 
Albrow M. D., et al., 1999, ApJ, 522, 1022 
Alcock C, et al., 1996, ApJ, 463, L67 
Alcock C, et al., 1997, ApJ, 479, 1 19 

Ansari R., 1994, Une methode de reconstruction photometrique 

pour r experience EROS, preprint LAL 94-10 
Beaulieu J.-P, et al., 2006, Nature, 439, 437 
Bennett D. P, Rhie S. H., 2002, ApJ, 574, 985 
Bensby T., et al., 2011, A&A, 533, A134 
Bessell M. S., Brett J. M., 1988, PASP, 100, 1 134 
Bond I. A., et al, 2001, MNRAS, 327, 868 
Bozza v., 2001, A&A, 374, 13 
Bozza v., 2010, MNRAS, 408, 2188 
Bramich D. M., 2008, MNRAS, 386, L77 
Bramich D. M., et al., 201 1, MNRAS, 413, 1275 
Burgdorf M. J., Bramich D. M., Dominik M., Bode M. F, Home 

K. D., Steele I. A., Rattenbury N., Tsapras Y, 2007, P&SS, 55, 

582 

Calchi Novati S., de Luca F, Jetzer P., Mancini L., Scarpetta G., 

2008, A&A, 480, 723 
Cassan A., 2008, A&A, 491, 587 
Cassan A., et al., 201 1, Nature accepted 
Chabrier G., 2003, PASP, 115, 763 

Choi J.-Y, et al., 2011, Characterization of Lenses and Lensed 
Stars from Extremely High-Magnification Events, to be submit- 
ted to ApJ 

Dominik M., 1995, A&AS, 109, 597 



Dominik M., 1998a, A&A, 333, L79 
Dominik M., 1998b, A&A, 333, 893 
Dominik M., 1999a, A&A, 341, 943 
Dominik M., 1999b, A&A, 349, 108 
Dominik M., 2006, MNRAS, 367, 669 

Dominik M., 2010, General Relativity and Gravitation, 42, 2075 

Dominik M., 201 1, MNRAS, 41 1, 2 

Dominik M., et al., 2002, P&SS, 50, 299 

Dominik M., et al., 2007, MNRAS, 380, 792 

Dominik M., et al., 2008a, in Y.-S. Sun, S. Ferraz-Mello, & 
J.-L. Zhou ed., lAU Symposium Vol. 249 of lAU Sympo- 
sium, ARTEMiS (Automated Robotic Terrestrial Exoplanet Mi- 
crolensing Search) - Hunting for planets of Earth mass and be- 
low, pp 35-41 

Dominik M., et al., 2008b, AN, 329, 248 

Dominik M., et al., 2010, AN, 331, 671 

Dominik M., Hirshfeld A. C, 1996, A&A, 313, 841 

Einstein A., 1936, Science, 84, 506 

Erdl H., Schneider P, 1993, A&A, 268, 453 

Fried D. L., 1978, JOSA, 68, 1651 

Gaudi B. S., 1998, ApJ, 506, 533 

Gaudi B. S., Han C, 2004, ApJ, 611, 528 

Gould A., 2004, ApJ, 606, 319 

Gould A., et al., 2010, ApJ, 720, 1073 

Gould A., Gaucherel C, 1997, ApJ, 477, 580 

Gould A., Loeb A., 1992, ApJ, 396, 104 

Grenacher L., Jetzer R, Strassle M., de Paolis F, 1999, A&A, 351, 
775 

Griest K., Hu W., 1992, ApJ, 397, 362 

Griest K., Safizadeh N., 1998, ApJ, 500, 37 

Grundahl F., Christensen-Dalsgaard J., Arentoft T., Frandsen S., 
Kjeldsen H., J0rgensen U. G., Kjaergaard P., 2009, Communica- 
tions in Asteroseismology, 158, 345 

Han C, 2007, ApJ, 661, 1202 

Han C, et al., 2009, ApJ, 705, 1116 

Han C, Gould A., 1995, ApJ, 449, 521 

Han C, Jeong Y, 1998, MNRAS, 301, 231 

Home K., Snodgrass C, Tsapras Y, 2009, MNRAS, 396, 2087 

Jaroszynski M., Paczynski B., 2002, AcA, 52, 361 

J0rgensen U. G., 2008, Physica Scripta, T130, 014008 

Kains N., et al., 2009, MNRAS, 395, 787 

Kervella P., Bersier D., Mourard D., Nardetto N., Fouque P., 

Coude du Foresto V., 2004, A&A, 428, 587 
Kiraga M., Paczynski B., 1994, ApJ, 430, LlOl 
Mao S., Di Stefano R., 1995, ApJ, 440, 22 
Mao S., Paczynski B., 1991, ApJ, 374, L37 
Muraki Y, et al., 201 1, ApJ, 741, 22 

Night C, Di Stefano R., Schwamb M., 2008, ApJ, 686, 785 

Paczynski B., 1986, ApJ, 304, 1 

Paczynski B., 1991, ApJ, 371, L63 

Paczynski B., 1996, ARA&A, 34, 419 

Sumi T., et al., 2010, ApJ, 710, 1641 

Tsapras Y, et al., 2009, AN, 330, 4 

Tubbs R. N., Baldwin J. E., Mackay C. D., Cox G. C, 2002, A&A, 

387, L21 
Udalski A., 2003, AcA, 53, 291 

Udalski A., Szymaiiski M., Kaluzny J., Kubiak M., Mateo M., 

1992, AcA, 42, 253 
Udalski A., Szymahski M., Kaluzny J., Kubiak M., Mateo M., 

Krzemihski W., Paczynski B., 1994, AcA, 44, 227 
Yoo J., et al., 2004, ApJ, 616, 1204 



© 0000 RAS, MNRAS 000, 000-000 



OGLE-2008-BLG-5 1 - weak microlensing anomaly 1 7 



(Affiliations continued) 

South African Large Telescope, P.O. Box 9 Observatory 7935, South Africa 

Department of Physics and Astrophysics, Faculty of Science, Nagoya University, Nagoya 464-8602, Japan 
^^Armagh Observatory, College Hill, Armagh, BT61 9DG, Northern Ireland, United Kingdom 

^'^ Korea Astronomy and Space Science Institute, 776 Daedukdae-ro, Yuseong-gu, Daejeon 305-348, Republic of Korea 
^^Istituto Intemazionale per gli Alti Studi Scientifici (IIASS), Vietri Sul Mare (SA), Italy 
'^^ Max Planck Institute for Astronomy, Kdnigstuhl 17, 6191 17 Heidelberg, Germany 
^'^INFN, Gruppo Collegato di Salerno, Sezione di Napoli, Italy 

Centra de Astro-Ingenierta, Departamento de Astronomta y Astrofi'sica, Pontificia Universidad Catdlica de Chile, Casilla 306, Santiago, Chile 

Deutsches SOFIA Institut, Universitdt Stuttgart, Pfaffenwaldring 31, 70569 Stuttgart, Germany 
^° SOFIA Science Center, NASA Ames Research Center, Mail Stop N2I1-3, Moffett Field CA 94035, United States of America 
^^Danmarks Tekniske Universitet, Institut for Rumforskning og -teknologi, Juliane Maries Vej 30, 2100 Kbenhavn, Denmark 

Bellatrix Astronomical Observatory, Via Madonna de Loco 47, 03023 Ceccano (FR), Italy 
^^Department of Physics, Sharif University of Technology, P. O. Box 1 1 155-9161, Tehran, Iran 
^^Institut d'Astrophysique et de Geophysique, Allee du 6 Aout 17, Sart Tilman, Bat. B5c, 4000 Liege, Belgium 
^^European Southern Obserx'citory (ESO) Alonso de Cordova 3107, Casilla 19001, Santiago 19, Chile 

Max Planck Institute for Solar System Research, Max-Planck-Str. 2, 37191 Katlenburg-Lindaii, Germany 
^''Astrophysics Group, Keele University, Staffordshire, ST5 5BG, United Kingdom 

^^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead CH4I ILD, United Kingdom 
^^Instituto de Astroflsica de Andalucfa, Glorieta de la Astronomia s/n, 18008 Granada, Spain 

'^^ Dark Cosmology Centre, Niels Bohr Institutet, K0benhavns Universitet, Juliane Maries Vej 30. 2100 K0benhavn 0, Denmark 

University of Canterbury, Department of Physics and Astronomy, Private Bag 4800, Christchurch 8020, New Zealand 
'^^ Department of Physics, 225 Nieiiwland Science Hall, University of Notre Dame, Notre Dame, IN 46556, Uiiited States of America 
'^^ McDonald Observatory, 16120 St Hwy Spur 78 #2, Fort Davis, TX 79734, United States of America 
"^"^ University of Tasmania, School of Mathematics and Physics, Private Bag 37, Hobart, TAS 7001, Australia 

'^^ Institute of Geophysics and Planetary Physics (IGPP). L-4I3, Lawrence Livermore National Laboratoiy, PO Box 808, Livermore, CA 94551, United States 
of America 

Department of Physics, Uiiiversity ofRijeka, Omladinska 14, 51000 Rijeka, Croatia 
^"^ Technische Universitdt Wien, Wiedner Hauptstr. 8-10, 1040 Wien, Austria 

'^^NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 South Wilson Avenue, Pasadena, CA 91125, United States of America 

Perth Observatory, Walnut Road, Bickley, Perth 6076, WA, Australia 

South African Astronomical Obserx'atory, P.O. Box 9 Observatory 7935, South Africa 
^^Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, United States of America 

Universidad de Concepcion, Departamento de Astronomia, Casilla 160— C, Concepcion, Chile 

Department of Physics, Texas A&M University, 4242 TAMU, College Station, TX 77843-4242, United States of America 

Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, United States of America 

Department of Physics, Ohio State University, 191 W. Woodruff, Columbus, OH 43210, United States of America 
^^Mt. John Observatory, P. O. Box 56, Lake Tekapo 8770, New Zealand 

School of Chemical and Physical Sciences, Victoria University, Wellington, New Zealand 
^^Department of Physics, Koncm Uiiiversity, Nishiokamoto 8-9-1, Kobe 658-8501. Japan 

Nagano National College of Technology, Nagano 381-8550, Japan 
^^Tokyo Metropolitan College of Aeronautics, Tokyo 116-8523, Japan 

Department of Earth and Space Science, Osaka University, l-I Machikaneyama-cho, Toyonaka, Osaka 560-0043, Japan 

Department of Earth Sciences, National Taiwan Normal University, No. 88, Section 4, Tlngzhou Road, Wemlmn District, Taipei 11677, Taiwan 



© 0000 RAS, MNRAS 000, 000-000 



